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ABSTRACT 


As part of a program to observe circumstellar emission from late type 
stars, lUE observations of R Aqr (M7 + pec) have been obtained in low 
dispersion. Strong permitted, semi -forbidden and forbidden emission lines 
are seen, superimposed on a bright ultraviolet continuum. We deduce that 
the strong emission line spectrum that involves C III], C IV, Si III], 

[0 II] and [0 III] probably arises from a dense compact nebula the size of 
which is comparable to the orbital radius of the binary system of which 
R Aqr is the primary star. The low excitation emission lines of Fe II, 

Mg II, 0 I and Si II probably arise in the chromosphere (T~ 10,000 K) of the 
R Aqr. The secondary is probably a white dwarf, comparable to or somewhat 
brighter than the sun, since such a star can produce enough ionizing photons 
to excite the continuum and emission line spectrum and yet be sufficiently 
faint as to escape detection by direct observation. We attribute the UV 
continuum to Balmer recombination from the dense nebula and not to blackbody 
emission from the hot companion. 

I. Data, Analysis, and Introduction 

Ultraviolet spectra (Figures 1 and 2) of the late type. star 
R Aqr (M7 + pec) have been obtained with the International Ultraviolet 
Explorer (lUE) and reveal intense emission lines and continuum. This is 
consistent with earlier optical observations of MerrillL* 3 which 
indicate that the system consists of hot stellar companion and a relatively 
cool late type star. In low dispersion the spectrum between 1200 A to 
3200 % shows strong permitted emission lines of C IV (1548 A, 1550 A), 

Si III (1883 A, 1892 A), C III (1907 8, 1909 A) and Mg II (2796 A, 

2803 R), forbidden emission lines of 0 II (2470 A) and probably [0 III] 
(2321 A). The UV observations are thus consistent with lines of [0 III] 
(4929 K, 5007 K) observed in the optical spectra of Merril r and [0 II] 

(3726 3729 K) of Ilovaisky and Spinrad^. 

The strong lines of He II, C II, C IV , 0 I, 0 III , 0 IV, S II, Si IV 
and Fe II are evident in the spectrum. The identification of N V and Si II 
1304 K, 1309 K is ambiguous because other lower excitation lines of nitrogen 
are not present, and similarly for silicon. Si II 1265 A is not observed. A 
number of Fe II features in various multiplets are also identified. 
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In addition to the bright emission line spectrum there is a general 
ultraviolet continuum, the intensity of which appears independent of wave- 
length over the spectral range observed. We attribute the origin of this 
continuum to hydrogen recombination rather than H I two-photon emission. 

Of particular interest in our observations is the distinct lack of a 
stellar UV continuum that should be present if an 0 or B type main sequence 
star ^ is the source of excitation in the nebula, as has been suggested by 
Merrill'^. We conclude that the continuum observed most likely originates 
not from a stellar source, but from a low excitatioocnebula with an electron 
temperature T'' 15,000 K and characteristic size 10 ^cm. The source 
excitation is a sub-luminous central planetary nebula star or bright white 
dwarf of T*> 50,000 K, whose orbit about the primary M7 star is comparable 
to the size of the ionized nebula, i.e. lOl*^ - 10^5 cm. The details of 
our conclusions are discussed in the following sections. 

The low excitation lines of Fe II and Mg II have been previously 
observed in the spectra of single late type stars by Carpenter and Wing 5, 
and the presence of the above lines as well as 0 I and Si II lines in the UV 
spectrum of R Aqr argues strongly for a cool chromosphere T*~10,000K for 
the primary M7 star . However, we assume from our analysis that 
other high excitation lines observea in our data do not arise from the 
companion directly. We attribute the formation of the majority of strong 
lines to a compact nebula that is excited by emission from the hot companion. 

The observed continuum flux can be used to obtain the general 
parameters of the nebula. These parameters can then be checked against 
those derived from our analysis of the continuum spectrum. For that purpose 
we have used the combined strengths of C II (2325, 2327, 2328 ^), C IV 
(1548, 1551 A), C III] (1907, 1909 A), [3 II] (2470 A), [0 III] (2321 A) and 
0 IV] (1400, 1401, 1405, 1407 A). We have selected these lines because they 
consist of various ion species, they are the strongest features in the 
spectrum, and because there is little ambiguity in identification. However, 
it is not clear if the 1402 A feature is due entirely to 0 IV] , and what 
portion of the broad 2328 A feature is due to [0 III] or C II. 

The observed continuum is essentially flat and rises slightly toward 
long wavelengths. The continuum, therefore, cannot be due to a star since a 
stellar continuum would vary with wavelength by more than an order of 
magnitude over the spectral range observed. Balmer recombination and H I 
two-photon emission arising from a nebula are possible mechanisms to explain 
this continuum. However, shortward of 3200 A the two-photon continuum-of 
hydrogen will dominate if the densities are sufficiently low (;S10^ cm"'^), 
and will produce a prominent peak in intensity around 1400 A (Bohlin, 
Harrington and Stecker^); this peak is not observed in our data, and we 
conclude that the two-photon process is not the dominant mechanism. On 
the other hand, the Balmer recombination continuum depends only weakly on 
temperature for T~15,000 K and the expected flux varies only by a factor of 
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two between 1200 to 3200 K, We conclude that this mechanism dominates and 
that the densities involved must be^lOS cm~3, since densities 
appreciably lower than this value would result in a conspicuous two-photon 
continuum peak around 1400 A. Given this lower limit in density we can 
calculate a corresponding upper limit for the diameter of an ionized nebula 
responsible for the observed flux, ~5 x 10^^ cm. 

A further argument in support of a nebula with this characteristic size 
and density as that found above can also be made from observations of 
emission lines in the visible by Merrill^. Merrill^ determined an 
orbital period for the hot companion of approximately 27 years from radial 
velocity measurements. The corresponding size of the semi-major axis for an 
elliptical orbit with this period is 1.7 x lO^^ to 2.1 x IQl^ cm, for a 
mass ratio of the primary and secondary of 1:1 and 3:1, respectively 
(assuming one solar mass for the secondary). If the ionized nebula was 
appreciably larger than the separation of the stars one would not expect to 
observe substantial variations in line strengths. However, Ilovaisky and 
Spinrad^ have compared their visual spectral data of R Aqr to earlier 
observations of Merrill and found no evidence for a hot stellar companion 
whatsoever. From this they suggest the emission properties of the spectrum 
are probably time-dependent. Merrill^ also found that the apparent 
position of the nebular emission appears to vary with time. Accordingly, 
this is consistent with our model that suggests the ionized nebula is 
comparable in scale to the size of the binary orbit. At a distance of 
260 pc a central ionized cloud of scale size L = 2 x lOl^ and electron 
density of ne = 1.5 x 107 cm-3 is sufficient to explain the observed 
recombination continuum. 

The densities we obtained above from observations of the UV continuum 
can also be roughly checked by the [0 III] line strengths observed by 
Merrill^ who found that the 4363 A line is unusually strong, and that 
prior to 1934 "the ratio of its intensity of 4959 A is equaled or exceeded 
in only one other nebula, IC 4997". The intensity of the hydrogen lines 
relative to those of [0 III] R Aqr is similar to that observed in most 
planetary nebulae, H3 being about equal to 4959 a prior to 1922. After this 
time the hot component started to dominate the spectrum for a few years 
reaching a maximum in 1933-34°. From Aller and tiller^ we 
estimate the intensity ratio of the 4363 A to 5007 A and 4959Ai.e. 

1(4363 A)/I(4959 K, 5007 A) 0.1 for the planetary nebula IC 4997 observed 
in 1922. Based upon the observations of IC 4997 (Aller and tiller^ and 
the statement of Merrill 3 concerning the similarity of this ratio of 
R Aqr to that of IC 4997 prior to 1934 (i.e. 1922), it follows that the 
strength of the 4363 A line agrees with densities lO^^ne^^lO^ cm"3 
for temperatures 3 x 10A<T^8 x 103, respectively (Kafatos and Lynchl^. 

Moreover, as deduced from our continuum measurements the nebular density 
cannot be more than ~103 cm"3 since the [0 II] and 0 [III] lines would 
be suppressed, and the 4363 A line would be even stronger. Additionally, if 
densities were much higher the nebula would be comparable in extent to the 
size of the primary M7 star, which is unrealistic on physical grounds. It 
follows that the nebular parameters from the foregoing arguments are in the 
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range 105<ne^l08 cm-3, 10l4<L<10l5 cm, and 104sT<3 x 
10^ K. We adopt a model ne< 107 cm-3, L~ 2 x lOl^ cm and 
T~ 15,000 K. 

b.) Emission Lines 

We have also deduced the nebular parameters and the relative ionic 
abundances by two different methods, i) by using the carbon line strengths, 
ii) by using the oxygen line strengths. We have assumed that the scale size 
of the nebula is L~ 2 x lOl^ cm; although the general trend of our results 
appears somewhat insensitive to this parameter, we have assumed normal 
cosmic abundances for the various elements under consideration!!. 

i) carbon line strengths (model A) 

The semi -forbidden lines of C II and C III and the allowed lines of 
Cp-IV can be used with one another to find the product n^ l3 '2, 13, 14, 

It is essential that we identify the 2328 A feature as C II], 
otherwise no self-consistent model can be constructed using the carbon lines 
(even if the [0 III] were present it would not be more than 0.1 of the 
total intensity of the feature). The results of this calculation are shown 
in Table 1 (Model A). The ionic abundance N(0 II) and N(0 III) can then be 
obtained (the latter is an upper limit since some 0 IV could also be 
present). We can also obtain the ionic abundance of He III from the 1640 % 
line using this carbon line strength analysis. However, the ionic abundance 
of the He III is found to be large in this case and all of the helium would 
have to be doubly ionized. We have deduced the H-a (6563 K) and [ 0 III] 

5007 A) as well as the flux of the continuum in Table 1 (Model A). Although 
Model A is not unique it does indicate a general trend in our data. 

Reasonable ionic abundances for the oxygen ions can be obtained if we use the 
carbon line strengths to n^ L^. However, the densities would have to be 
generally lower than might be suggested by the strong [0 III] 4363 a line 
observed by Merrill 3. Additionally, the computed flux level of the 
continuum from Model A is too low to explain our data. 

ii) oxygen line strengths 

Assuming the 2328 A feature is mainly due to [0 III] and the 1402 ^ 
feature is mainly 0 IV], one can obtain values of ne and t!^. Since 
collisional depopulation becomes important above 106 cm"3, the 
forbidden lines values of n^ obtained from a single temperature are not 
unique, but are relatively insensitive to temperature. A typical case is 
shown in Table 1 for T = 15,000 K and L = 2 x lo!'^ cm (Case B). The 
continuum deduced from the line strengths of oxygen agree well with the 
observed continuum. Moreover, the deduced ionic abundance of He III agrees 
with the ionic abundance of 0 III (helium is essentially singly ionized, 
whereas oxygen is mostly doubly ionized). 

On the other hand, if we use the cosmic abundances of carbon in this 
analysis, the line strengths should be a factor of ~50 larger than what is 
observed. The only alternative is to assume that atomic carbon is under- 
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abundant by a factor ^ 50 in the ionized nebula. The depletion of carbon 
could be the result of the precipitation of carbon into grains. The 
relatively low abundance deduced in our analysis appears generally to be 
the case since this result will not change even if parameters in Table 1 are 
varied (say the temperature is varied- between 10,000 K to 15,000 K). At 
this point, however, it is not possible to distinguish between either Models 
A or B. It suffices here to say that Model B appears more attractive since 
it does account for the continuum. 

It is of interest to note that the values of the nebular parameters 
deduced in our analysis, i.e. L = 2 x loi^ cm, ne~10^ - 10^ cm“3 
and T~15,000 K agree with the general parameters for nebular emission in 
symbiotic stars of^®. 

The compact nebula could be entirely the result of mass loss from the 
primary star. Applying the equation of continuity and estimating the escape 
speed of the M7 giant to be 24 km s"S that was obtained using the 
period-density relation for a period = 387 days and an assumed stellar mass 
of.lMQ to 3Mo , we find M = 10“7 M^yr-l for a nebula of radius rEL/2 
lo'^ cm and density ng~1.3 x 10' cm”^. This mass loss rate is 
probably a lower limit since, as will be shown later on, the hot companion 
is too faint to ionize the entire nebula. It is also unlikely that all the 
material lost by the star is still ionized. 

c.) properties of the stellar companion 

We have already seen that the observed continuum cannot be due to a star 
and most likely arises from a nebula. However, since we have assumed that 
the source of nebular excitation is a sub-luminous white dwarf, the nebular 
continuum flux observed of 10“13 ergs cm"2 s“^ A"1 places an upper 
limit of flux that is contributed by the companion. In Table 2 we show the 

stellar parameters for the unseen companion if its flux contribution in the 

continuum is this upper limit. 

In the first column of Table 2 we assume a temperature for the hot com- 
panion. The second column gives the corresponding stellar radius ifthe 
continuum at the detector is 10"l3 ergs cm"^ s"t p^-1 at 1200 A. The 
third column gives the luminosity in solar units. Columns 4 and 6 give the 
absolute flux and apparent visual magnitudes, respectively. The last column 
indicates the number of ionizing photons Ni(s~^) emitted by the star. 

This upper limit of continuum flux suggests a star whose apparent 
magnitude is too faint to be observed today. However, in 1934 it attained 
mv~8 magnitudes. In order for it to be observable today the continuum 

would have to be lO'^ times greater. The 1934 event appears, therefore, to 

have been en eruption of the hot companion that was not sufficiently strong 
to be classified a nova. It is possible that this eruption was triggered by 
mass transfer from the primary to the secondary. 


359 



We also note that the companion can ionize the dense nebula in the 
system. The stellar temperatures required are in the approximate range of 
10» to 1.5 X 105 K, although temperatures as low as 5 x 10^ K could be 
assumed if the density of the nebula was slightly lower than that shown in 
Model B. It is most likely the case that the compact ionized nebula is 
"ionization" bound rather than "density" bound and therefore, the mass rate 
for the primary of 10-' M© yr-1 is likely a lower limit. We also find 
that nebular temperatures greater than 15,000 K are required, otherwise the 
density implied by the analysis based on the oxygen line strengths would be 
high, and a much brighter star would be required to ionize the nebula. 

We find that the nebular parameters deduced from our oxygen line 
analysis, ng^lO^ cm"3, 1.5 x 104<T<2.5 x 10l4 cm, and the companion 
5 x 10^^T;^1.5 x 105 K, 0.7 7L can account for the nebular 

continuum, explain the observed emission lines, and provide sufficiently 
low stellar luminosity that the companion could not be seen directly. 

III. Summary 

The UV emission observed in our data and the forbidden line emission 
observed by Merrill^ and Ilovaisky and Spinrad^ are most likely the 
result of the excitation of a nebula by a white dwarf. If an 0 or B type 
main sequence star is postulated as the excitation source of the nebula, 

such a model could not reconcile the continuum properties of the ultraviolet 

spectrum. Accordingly, our model is constrained to adopt a model in which a 
bright white dwarf (few Lq) is the companion to the M7 giant. The fact 

that a white dwarf is capable of supplying enough ionizing photons to excite 

the emission lines observed further strengthens this interpretation. We can 
summarize the general properties of our model as follows: R Aqr is a 
symbiotic star system that most likely consists of an M7 primary and white 
dwarf companion. The 27 year period we have adopted from Merrill^ for the 
companion star is such that for reasonable mass ratios of 1:1 and 1:2 
(assuminq the dwarf to be a IM© ) the physical separation of the stars is a 
few X lOl^ cm, which is also the approximate dimensions of the ionized 
nebula. The faint hot companion star is itself not sufficiently luminous to 
be observable directly. Its presence, however, is manifested in the ionizing 
effects which it has on its immediate surroundings, which create a low 
excitation nebula. Further observations in the radio, visible and ultra- 
violet would be useful in monitoring the time- dependence of the different 
emitting regions. 
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Figure 1: 

Low dispersion (~6 A resolution) spectrum of R Aqr obtained on 

July 25, 1979 in the wavelength range 1200 A to 2000 R using the large, 
aperture (10" x 20") of the IDE spectrometer. The exposure time on the 
SEC Vidicon camera was 10 minutes. The Lyman-a line at 1216 R has a 
combined intensity of the geocoronal and stellar 
Lyman-a emission. The spectrum is deconvolved and geometrically and 
photoelectrically corrected. The absolute flux scale is accurate to a 
factor two using the standard lUE reduction program. 
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F i QU r*G 2 • 

Low dispersion spectrum of R Aqr between 2000 K and 3200 A obtained with the 
large aperture on lUE. Exposure time was 20 minutes. The feature 
identified with a question mark is probably an artifact in the detector. 

The continuum level is seen to rise slightly toward long wavelengths. 
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TABLE 1. 


Nebula P^ameters 
L = 2 X 10^^ cm, T = 15.000K 


Case A 

Carbon Line Analysis 


, 3, -3. 
n^ L (cm ) 

4.4 X 10^® 

N(C II) 

0.53 

N(C III) 

0.32 

N(C IV) 

0.15 

Hg (cm'^) 

2.3 X 10® 

N(0 II) 

0.81 

N(0 III) 

0.19 

N(He III) 

1 

Flux (H-oC)_-j 6563A 
(ergs cm”2 s' ) 

1.2 X 10'^^ 

Flux Co III] 5007A 
(ergs cm"2 s"1 ) 

1.5 X 10'^^ 


10-1® 


(ergs cm"*^ s' A' ) 


Case B 

Oxygen lire Analysis 


Hg (cm'^) 

1.25 X 10®^ 

N(0 II) 

0.06 

N(0 III) 

0.86 

N(0 IV) 

0.08 

ng(cm'^) 

o 

r— 

>< 

tr> 

CVi 

N(C II)* 

0.09 

N(C iin* 

0.55 

N(C IV)* 

0.36 

N(He III) 

0.02 

Flux (H-oC) 
(ergs cm"2 s~^) 

3.3 X lO'll 

Flux [0 III ] 5007A 
(ergs cm"2 s‘>) 

5.5 X lO'll 

'^Contfnuun,<2'’°“> 
(ergs cm"2 s" ) 

2.0 X lO'l^ 


★ 


see discussion in text 
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TABLE 2 

HOT COMPANION 


T*(K) R*/R© L^/Lg> My my 


1.5 X 10^ 4.1 X 10“^ 7.4 

10^ 5.6x10“^ 2.8 

5.0x10^ 1.2x10"^ 0.75 

Tj, : stellar temperature 

Rj. : stellar radius 

: stellar luminosity 

M^ : absoulte visual magnitude 

m^ : apparent visual magnitude 

N. : number of ionizing photons 


12 

19.1 

4.3 

X 

10*'* 

11.4 

18.5 

2.0 

X 


10.1 

17.2 

5.0 

X 

10« 


emitted by companion 



